Multi-wavelength spectroscopy of the bipolar outflow from 

CepheusE ^ 
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ABSTRACT 

Cepheus E is the site of an exceptional example of a protostellar outflow with a very young 
dynamical age and extremely high near infrared luminosity. We combine molecular spectroscopic 
data from the submillimeter to the near infrared in order to interpret the rotational excitation of 
CO and the ro- vibrational excitation of H2. We conclude that C-type shocks with a paraboloidal 
bow shock geometry can simultaneously explain all the molecular excitations. Extinction accounts 
for the deviation of the column densities from local thermodynamic equilibrium. A difference in 
the extinction between the red and blue-shifted outflow lobes may account for the measured flux 
difference. The outflow is deeply embedded in a clump of density 10^ cm~^, yet a good fraction 
of atomic hydrogen, about 40%, is required to explain the excitation and statistical equilibrium. 
We propose that this atomic component arises, self-consistently, from the dissociated gas at the 
apex of the leading bow shocks and the relatively long molecule reformation time. At least 20 
bow shocks are required in each lobe, although these may be sub-divided into smaller bows and 
turbulent shocked regions. The total outflow mechanical power and cooling amounts to over 
30 Lq, almost half the source's bolometric luminosity. Nevertheless, only about 6% of the clump 
mass has been set in outward motion by the outflow, allowing a collapse to continue. 



Subject headings: Shock waves - Molecular processes - ISM: jets and outflows 
dynamics - ISM: molecules - stars: mass-loss 



ISM: kinematics and 



1. Introduction 

Twin coUimated lobes of molecular gas stream 
away from newly forming stars (Bachiller 1996). 
These bipolar outflows possess particularly high 
power and thrust during the main phase of inflow 
onto the protostar. The driving agents are often 
recognised as pulsating supersonic jets, originat- 
ing from near the protostellar surface. The ex- 
tended environment is pushed and shocked, pro- 
ducing bow-shaped structures called Herbig-Haro 



^ Based on observations with ISO, an ESA project with 
instruments funded by ESA Member States (especially the 
PI countries: France, Germany, the Netherlands and the 
United Kingdom) and with the participation of ISAS and 
NASA. 



(HH) objects (Rcipurth & Bally 2001). The HH 
objects are often optically invisible because the 
protostar is deeply enshrouded in a dusty cloud. 
We analyse here one such outflow, originating 
from the Class source Cepheus E - MM (Lefloch 
et al. 1996), through combined near-infrared data, 
mid- and far-infrared (FIR) ISO spectra and sub- 
millimeter maps. The numerous emission line 
strengths at these long wavelengths constrain and 
relate the gas components. 

CepheusE contains a powerful outflow from a 
luminous source. A CO-derived kinetic luminos- 
ity of 0.2 Lq (Moro-Martfn et al. 2001), a far in- 
frared line luminosity of 2.8 Lq (Giannini et al. 
2001) and an H2 ro-vibrational line luminosity of 
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0.7 L0 (Frocbrich ct al. 2003) have been derived, 
based on simple assumptions and an estimated 
distance of ~730pc. The driving protostar has 
a bolomctric himinosity of ^ 80 Lq (Lcfloch ct al. 
1996; Froebrich et al. 2003) and is surrounded by a 
protostellar envelope of 25 Mq (Lefloch et al. 1996; 
Chini ct al. 2001). This suggests that the driving 
young stellar object will develop into an interme- 
diate mass object, consistent with a straightfor- 
ward evolutionary model (Frocbrich ct al. 2003). 
This contrasts with other well-known text-book 
examples of jet-driven Class outflows, such as 
HH211 and HH212, which may be powered by 
low-mass or solar-mass stars. These statements 
are, however, based on evolutionary assumptions 
which need to be tested. 

The outflow is apparently driven by jets or bul- 
lets with radial velocities of -120 and -|-80kms~^ 
(Smith, Suttner, & Yorkc 1997), as derived from 
CO spectra (Lefloch et al. 1996; Hatchell et al. 
1999). The proper motion of an optical knot at 
the edge of the southern blue-shifted lobe, HH 377, 
is 107±14kms~^ (Noriega-Crespo and Garnavich 
2001) with a radial velocity of -70±10kms-^ (Ay- 
ala et al. 2000). This implies that the jet is 
more dense than the environment through which 
it propagates. Reversing the well-known formula 
for thrust balance yields a jet-ambient density ra- 
tio of ?7 = l/(vjet/vfeouj — 1)^ = 2.0 (taking the 
radial velocity components and assuming the gas 
ahead of the bow shock is stationary). The inclina- 
tion angle to the line-of-sight is ~tan~^ (107/70) 
= 57°±7° and the bow speed is 128±12kms-^ 

We feature CepE in this study because of its 
extreme youth and high luminosity. Several other 
properties have attracted attention to this outflow. 

• The low excitation. Both near-infrared H2 
and optical H lines yield remarkably low 

excitation. Despite the strong emission 
fluxes stemming from vibrationally excited 
H2 states, the excitation temperature is low, 
with its precise value sensitive to the exci- 
tation levels between which it is measured 
(Eisloffel et al. 1996). The shock speed of 
HH 377 derived from H/3 and Ha fluxes are 
under 20 km s^^ (Ayala et al. 2000). 

• The constant molecular excitation. The 
H2 line ratios vary remarkably little over 
the entire outflow on small and large scales 



(Eisloffel et al. 1996). 

• The evidence for a precessing underlying 
flow (Eisloffel et al. 1996) 

We here analyse the radiative shock waves, si- 
multaneously modelling the many strong emis- 
sion lines from carbon monoxide and hydro- 
gen molecules through both rotational and ro- 
vibrational transitions. To accomplish this we 
combine ISO and ground-based data. We thus 
hope to go a stage further than previous stud- 
ies of K-band spectra and ISO data, analysed 
in isolation. Secondly, our model includes bow 
shock dynamics and ambipolar diffusion (C-shock) 
physics, both of which have not yet been applied 
to recent Cepheus E spectroscopic data (although 
Ladd & Hodapp (1997) compared Smith's (1995) 
tabulated C and J-type planar models to K-band 
spectra). 

Using refurbished shock codes, tested by Froe- 
brich et al. (2002a) for the outflows from CepA 
and L 1448, we model the outflow as planar and 
bow shocks with a symmetric shape Z oc (in 
cylindrical coordinates). The main aim of this 
process is to draw conclusions about the physics, 
chemistry and dynamics of the shocks and the 
properties of the surrounding gas. This analysis 
will also provide hints concerning the early evolu- 
tion of an outflow and possibly of the source itself. 

In Sect. 2 we describe the observations and data 
reduction. The main results are presented in 
Sect. 3. The results of the extinction modelling 
(Sect. 4), and modelling of the H2 and CO line 
fluxes (Sect. 5) are then given. 

2. Observations and Data Reduction 
2.1. Near infrared data 

The NIR spectra were obtained in the period 
26-29 August 1996 with the UH 2.2-m telescope at 
the Mauna Kea Observatory. The KSPEC spec- 
trograph (Hodapp et al. 1994), a cross-dispersed 
Echelle spectrograph, was employed. It is designed 
to provide medium-resolution spectra in the 1- 
2.5 Aim region. A HAWAII 1024x1024 detector 
array was used. We observed at three different 
positions in the outflow, as shown in Fig. 1. Data 
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Fig. 1. — Positions of the apertures for LWS ob- 
servations (circles) in Cepheus E superimposed on 
an image in the 1-0 S(l) line of molecular hydro- 
gen at 2.122 /xm. The slit positions of the three 
KSPEC spectra are also indicated. 

reduction, including flat-fielding, sky subtraction 

and extraction of the spectra, was done using our 
own MIDAS routines. Absolute flux calibration 
was not possible due to non-photometric weather 
conditions. The relative fluxes at wavelengths 
larger than A^2.4/Ltm (1-0Q() lines) are unre- 
liable because the flat-field was quite poor. H2 
emission lines were detected in the H and K band, 
which thus required wavelength calibration. For 
this purpose we used the OH night-sky emission 
lines and the tables of Rousselot et al. (2000). 

2.2. ISO data 

We employed three sets of ISO satellite data 
(Kessler et al. 1996). CepE was observed at 
two different positions with the Long Wavelength 
Spectrometer (Cesarsky et al. 1996) in the LWS 01 
grating mode. The apertures are indicated as cir- 
cles in Fig. 1 In the LWS 01 mode, a spectrum from 
43 to 197 /xm was obtained with a spectral resolu- 
tion of about 200. (See the ISO Handbook, Vol- 



Table 1: Observation log for Cep. E ISO data. 

TDT Object a (J2000) 5 (J2000) AOT 

56600912 CopES 23 03 13 +61 41 56 LWSOl 
56601113 Cep EN 23 03 13 +61 42 59 LWSOl 
79200740 CepE 23 03 13 +61 42 27 CAM04 



ume IV: LWS — The Long Wavelength Spectrom- 
eter^ and Cesarsky et al. (1996) for instrument 

and Astronomical Observing Template details). 

In addition, the region was imaged with the 
ISOCAM instrument (Cesarsky et al. 1996) 
through the CVF-filters in a wavelength range 
from 5 to 17 /im and a pixel size of six arcseconds. 
We refer to the work of Moro-Martfn et al. (2001), 
where these data were first presented. The ob- 
served field of the observations is larger than the 
NIR image shown in Fig. 1. (See the ISO Hand- 
book, Volume HI: CAM - The ISO Camera^ and 
Cesarsky et al. (1996) for instrument and AOT 
details.) A log of the ISO observations used here 
is given in Table 1. 

Data reduction was performed using ISO soft- 
ware (ISAP 1.6a and LIA 7.3 for the LWS data 
and CIA for the ISOCAM-CVF data), and the 
data from standard pipeline 8.7 (LWS) and 10 
(ISOCAM-CVF). In the LWS spectra glitches due 
to cosmic ray hits were removed as well as the 
heavy fringing which occurred in the spectra. Flux 
measurements were extracted using Gaussian fits 
to the lines plus second order polynomials for the 
background. Lines with a FWHM significantly 
smaller than the instrumental profile are not con- 
sidered as real, and are excluded from further anal- 
ysis. 

2.3. Submillimeter observations 

Submillimeter data of the CepE outflow were 
obtained in the ^^CO(3-2) rotational line. The 
observations were carried out in the period 18 - 21 
June, 1997 on the JCMT at Manua Kea. We em- 
ployed the B3 and C heterodyne receiver. Data re- 
duction was done using SPECX software. A spec- 
trum was measured at each point on a 7x11 grid 
with a grid point separation of six arcseconds. We 

■^http:/ /www.iso.vilspa.csa.cs/manuals/HANDBOOK/IV/lwsJib/ 
^http;/ /www. iso.vilspa.esa.es/manuals/HANDBOOK/III/cam_hb/ 
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also obtained some data for the lines ^^CO(4-3), 
CS(7-6) and SiO(8-7) which we briefly describe 
below. 

3. Results 

3.1. Near infrared data 

Almost 30 ro- vibrational transitions of molecu- 
lar hydrogen were detected in the H and K bands. 
These include the 1 S branch up to the 1- 
0S(9), the 1-0 Q branch up to 1 0Q(8), the 
2 - 1 S branch and even a few 3 2 S branch lines. 
Some lines were not detected because they are sit- 
uated in atmospheric absorption bands. Table 2 
provides the complete list of all relative fluxes de- 
tected with KSPEC. The line ratios are quite sim- 
ilar in all three positions. As we will find, the K- 
band extinction is quite high and might also have 
significant influence on the measured excitation of 
the molecular hydrogen. 

There are now 4 independent determinations 
of the flux ratio R21 = 1(2-1 S(l))/I(l-0 S(l)), 
at various locations. For location Na, we find 
R21 = 0.12, in agreement with Ayala et al. (2000) 
(their North location). We estimate R21 ~ 0.1 
from Ladd & Hodapp (1997), while Eisloffcl ct al. 
(1996) found R21 = 0.078±0.005 at location ND 
(roughly a 2" box). Hence the very low excitation 
reported in the latter work is not confirmed for 
this or the other locations. For the location in the 
southern lobe, SI, we find R21 = 0.13, while Ayala 
et al. (2000) find R21 = 0.11. This is consistent 
with typical 10% errors in the 2-1 measurements. 

A subset of the NIR emission lines detected here 
were also measured by Ayala ct al. (2000) for two 
positions, Na and SI, close to ours. This pro- 
vides an indirect means of calibrating our data. 
The main problem, however, is to determine the 
total flux of the K-band lines within the ISO 
apertures in order to directly compare and model 
the emission levels. The integrated 1 S(l) 
lobe fluxes are 1.4 x 10~^^ W m~^ (North) and 
2.7 X 10-15 W m-2 (South) (Eisloffel et al. 1996). 
The flux ratio is reasonably close to the value of 
0.44 measured at the location of the peak fluxes 
by Ayala et al. (2000). 



Table 2: Relative fluxes of H2 lines in the H- and K- 
band measured with KSPEC in the three observed 
slits (see Fig. 1). The fluxes are normalised to the 
1 - S(l) line. The errors are about 10 % and 25 % 
(for A^2.4/im). 



Transition 


An 


CepE 


CepE 


CepE 






[//ml 


Na 


Nd 


SI 


1 


-0S(9) 


1.6873 


0.025 


0.014 


0.035 


1 


-0S(8) 


1.7143 


< 0.015 


0.012 


0.026 


1 


-0S(7) 


1.7475 


0.097 


0.093 


0.144 


1 


-0S(6) 


1.7876 


0.063 


0.061 


0.090 


1 


-0S(5) 


1.8353 


0.261 


0.293 


0.322 


1 


-0S(4) 


1.8914 


0.099 


0.109 


0.351 


1 


-0S{2) 


2.0332 


0.328 


0.323 


0.349 


1 


-0S(1) 


2.1213 


1.000 


1.000 


1.000 


1 


-0S(0) 


2.2226 


0.339 


0.309 


0.281 


1 


-0Q(1) 


2.4059 


1.101 


1.010 


1.106 


1 


-0Q{2) 


2.4128 


0.417 


0.379 


0.341 


1 


-0Q{3) 


2.4231 


1.301 


1.160 


1.215 


1 


-0Q(4) 


2.4368 


0.367 


0.323 


0.340 


1 


-0Q(5) 


2.4541 


0.584 


0.648 


0.779 


1 


-0Q(6) 


2.4749 


0.276 


0.164 


0.209 


1 


-0Q(7) 


2.4993 


0.548 


0.4.30 


0.464 


1 


-0Q(8) 


2.5270 


< 0.044 


0.081 


0.136 


2 


-1S{7) 


1.8523 


< 0.033 


< 0.010 


0.051 


2 


-1S(4) 


2.0035 


< 0.016 


0.015 


0.040 


2 


-1S(3) 


2.0729 


0.084 


0.085 


0.095 


2 


-1S(2) 


2.1536 


0.043 


0.041 


0.048 


2 


-1S(1) 


2.2471 


0.108 


0.103 


0.117 


2 


-1S(0) 


2.3550 


0.019 


0.016 


0.027 


3 


-2S(6) 


2.0130 


< 0.011 


< 0.010 


< 0.007 


3 


-2S{5) 


2.0650 


0.016 


< 0.010 


0.015 


3 


- 2 S(4) 


2.1274 


< 0.016 


< 0.010 


< 0.008 


3 


- 2 S(3) 


2.2008 


< 0.010 


0.023 


0.025 


3 


-2S(2) 


2.2864 


< 0.010 


0.009 


0.015 


3 


-2S(1) 


2.3858 


< 0.010 


0.013 


0.027 



3.2. ISO LWS data 

The ISO-LWS spectra of the two outflow lobes 

show a wide variety of atomic fine structure and 
molecular lines. We detected rotational CO tran- 
sitions from Jup = 14-21 and some water lines in 
both lobes. The CO emission is on average 1.5 
times stronger in the southern outflow lobe, sug- 
gesting a higher filling factor of the beam or a 
higher CO abundance. The same applies for the 
H2O lines. All detected lines and fluxes in the 
LWS spectra are listed in Table 3. The LWS spec- 
tra for both outflow lobes can be found in Fig. 2 of 
Moro-Martfn et al. (2001). We flnd no particular 
anomalies with their derived fluxes. 
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Table 3: Observed lines in CepheusE North and 
South. The H2 fluxes measured in the ISOCAM 
data are co-added in the LWS apertures. The LWS 
fluxes are in IQ-^*^ Wm'^, the ISOCAM H2 fluxes 
are in 10"^^ Wm-^. 



EIgius nt 


Tr cHisiijion 


^0 


Ccp E 


Ccp E 






llJ,m] 


IN omi 


Q^ii+ii 


H2 


0-0 S(7) 


5.510 


4.2±0.8 


3.4±0.7 


H2 


0-0 S(6) 


6.107 


2.2±0.4 


2.0±0.4 


H2 


0-0 S(5) 


6.908 


8.3±1.7 


5.3±1.1 


H2 


0-0 S(4) 


8.023 


2.8±0.6 


3.6±0.7 


H2 


0-0 S (3) 


9.662 


6.1±1.2 


6.6±1.3 


H2 


0-0 S (2) 


12.275 


4.6±0.9 


5.8±1.2 


[OI] 


3Pl-3P2 


63.184 


93.4±5.0 


129±6.5 


0-H20 


321 -2l2 


75.. 380 


17.7±3.4 


<10.0 


0-H20 


5o5 -4i4 


99.492 


11.2±1.0 


<10.0 


CO 


25-24 


104.445 




<10.0 


0-H20 


221 - lio 


108.073 


)<15.0 


)<16.0 


CO 


24-23 


108.763 


CO 


23-22 


113.458 


}<22.4 


}<12.9 


0-H20 


4l4 -3o3 


113.537 


CO 


22-21 


118.581 


<10.0 


<13.2 


CO 


21-20 


124.193 


4.4±2.0 


<13.4 


P-H20 


4o4 -3l3 


125.353 


<8.0 




P-H20 


331 -322 


126.713 


<6.0 


<6.5 


CO 


20-19 


130.369 


6.5±1.1 


10.3±1.8 


CO 


19-18 


137.196 


9.2±1.0 


13.3±2.5 


P-H20 


3l3 -2o2 


138.527 


<2.5 


5.9±3.3 


CO 


18-17 


144.784 


12.3±1.5 


16.8±1.6 


[OI] 


3Po-3Pl 


145.525 


<2.0 


<3.0 


CO 


17-16 


153.267 


11.9±2.2 


19.8±1.0 


[CII] 


^P3/2-^Pl/2 


157.741 


72.4±1.1 


86.1±2.0 


CO 


16-15 


162.812 


15.7±1.2 


20.7±2.2 


CO 


15-14 


173.631 


19.9±3.0 


30.9±5.8 


0-H20 


3o3 -2i2 


174.626 


8.3±2.1 


14.0±0.9 


0-H20 


2l2 - loi 


179.527 


26.3±2.6 


29.1±2.2 


CO 


14-13 


185.999 


23.8±3.2 


28.9±3.5 



3.3. ISOCAM CVF data 

Pure rotational lines of the ground vibrational 
level of H2 were detected (0-0S(2) .. S(7)) with 
ISOCAM. Images in these lines as well as spectra 
of selected pixels are presented in Moro-Martm et 
al. (2001). To achieve our aim, to simultaneously 
model the CO and H2 lines, we co-add the spectra 
of all pixels of the ICOCAM images which are sit- 
uated in the LWS beam. Thus wc compare fluxes 
measured in the same aperture for both species, 
CO and H2. The co-added fluxes of the 0-0 S 
lines in the two LWS apertures are given in Ta- 
ble 3. 




R.A. offset(arcsec) from 23 01 10.£ 



Fig. 2. — Gray-scale and contour map of the in- 
tegrated flux in the CO 3-2 line of the CepheusE 
outflow. 



3.4. Submillimeter data 

The two outflow lobes are also present in the 
C0(3 2) map (Fig. 2). Note that the two lobes 
are elongated in a direction mis-aligned from the 
overall outflow direction. At the brightest points, 
the line shape shows signs of being optically thick 
(self absorption and rc-cmission). High velocity 
bullets with Vrei «100kms~^ are also detected. 
In the CS(7- 6) map we found weak emission only 
near the source position and extremely weak emis- 
sion at the three brightest knots of the outflow 
(SI, Na and Nd). Nevertheless, the outflow struc- 
ture is visible in this line. At the three positions 
where the CO (4 -3) line is observed we see the 
same line structure as for the CO(3-2) line: op- 
tically thick in the line centre and very fast CO 
bullets (vre; «100kms~^). The northern posi- 
tions observed in SiO(8-7) show no sign of a line. 
Only at the two positions in the south very weak 
emission is detected. 
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4. Extinction 
4.1. Methods 

To determine the luminosity and excitation of 
a deeply-embedded outflow it is necessary to re- 
move the extinction. Extinction reduces observed 
fluxes from their intrinsic values, especially for 
those from the shorter wavelength KSPEC data. 
Hence, in order to test shock models, we first ad- 
just the KSPEC fluxes. 

Since the CepE outflow is embedded in the 
parental cloud of the source (only the southern- 
most part is visible at optical wavelengths as 
HH377), the extinction is high. Moreover, the ex- 
tinction may vary not only over the field of the 
LWS aperture but also through a lobe. That is, 
the K-band data will be more representative of 
emission arising from low-excitation regions while 
ISO and CO data will sample more evenly. Here, 
however, since we do not have an extinction map 
and cannot distinguish components of contrasting 
extinction, we determine the extinction for each 
KSPEC location and then apply one value per 
LWS and ISOCAM beam. 

Three means are at our disposal to estimate 
the extinction. First, 1.25 mm contimuun emis- 
sion from cool dust is spatially coincident with 
the outflow (Lefloch et al. 1996). This can be in- 
terpreted as an H column of 6 x 10^^ cm in the 
lobe locations, corresponding to 3-4 mag of K- 
band extinction (incorrectly converted by Lefloch 
et al. (1996)). If half of this emission lies on av- 
erage in front of the lobe, an extinction of 1.5- 
2 mag would result. Submillimeter continuum ob- 
servations yield a mean H2 column of 3 x lO^^cm 
within a 30" radius, consistent with the above re- 
sult (Chini et al. 2001). 

Second, the Q-branch 1-0 lines beyond 2.4/xm 
arise from the same upper energy levels as the S- 
branch 1-0 lines within the K-band. Hence, differ- 
ential extinction between these wavelengths can, 
in theory, be determined exactly since the number 
of emitted photons must be proportional to their 
Einstein coefficients (radiative decay rates). Ab- 
solute extinction can then be derived from the dif- 
ferential value provided the properties of the dust 
are known. 

To apply this second method, we adopt a 
differential extinction of the form dlog{Ax) = 



0.4^K((2.12^m)/A)i '^ - 1.0). Using just the 
1-0 S(l) and 1-0 Q(3) lines (both originating 
from the same upper energy level), we determine 
K-band extinctions of Ak ^ 2.6 3.2 mag for 
the three KSPEC slit positions. In contrast, the 
Q(2)/S(0) ratio yields a uniform extinction of just 
l.libO.l mag. 

The Q-branch lines are, however, greatly ef- 
fected by transmission and are often inaccurate. 
One good test for accuracy is to inspect the ra- 
tio of the Q(3)/Q(l) lines which should be within 
the range [0.91,1.12] if filter or transmission dis- 
tortions are low and if the gas is excited within 
the range [1500K,3000K] (Smith 1995). This was 
not case for the previously derived values in Ladd 
& Hodapp (1997) and Moro-Martfn et al. (2001). 
Moro-Martfn et al. (2001) find values of 0.58 and 
0.86 for locations in the two lobes and Ladd & 
Hodapp (1997) find 0.78 and 0.88. Differential ex- 
tinction would only exacerbates the inconsistency. 
Here, we derive Q-branch Q(3)/Q(l) ratios of 1.18, 
1.16 and 1.10 for the 3 slits, which are comfort- 
ably consistent when some differential extinction 
is taken into account. Nevertheless, even in this 
case, extinctions derived from the S(l)/Q(3) and 
S(0)/Q(2) ratios do not agree. 

The third method also employs differential ex- 
tinction but makes use of a large fraction of the 
data. We plot the H2 columns from all the lev- 
els across the H and K bands and then determine 
the extinction range which best correlates the ac- 
cumulated data. To achieve this result, however, 
one must apply the normalised Column Density 
Ratio (CDR) method. Plotting the derived ab- 
solute columns of H2; Nj against the upper en- 
ergy level, Ej hides the information in the data 
since columns are spread out over three orders of 
magnitTide while error bars are only 10% on many 
data points. In the CDR method, the columns 
are normalised to the 1-0 S(l) fiux and divided 
by the equivalent LTE column of a gas at 2000K 
(see Froebrich et al. (2002a) and Eisloffel, Smith, 
& Davis (2000) for details). The results are shown 
in Fig. 3. 

Fig. 3 confirms that high extinction is present. 
Extinction is responsible for the distribution of the 
data points along the edges of a rhombus (top pan- 
els), which collapses to a line on applying the ex- 
tinctions indicated. The four sides of the rhom- 
bus consist of (i) the K-band 1-0 S-branch, (ii) 
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Fig. 3. — Column density ratios for the near infrared H2 ro-vibrational data, before applying extinction (top 
panels) and after adjustment for extinction (lower panels). The slit locations and extinction are indicated on 
each panel. Upper bounds are displayed as vertical lines. Squares represent 1-0 lines, with lighter squares 
denoting the less-reliable columns derived from the Q-branch, diamonds indicated 2-1, and triangles indicate 
CDRs from 3-2 lines. 



1-0 Q-branch, (iii) 2-1 S-branch and (iv) H-band 
1-0 S-branch. We note no evidence for UV excita- 
tion and fluorescent emission (as usually identifi- 
able from rotational transitions within the higher 
vibrational levels). Hence, the CDRs from the 2-1 
lines should lie level or below those from the 1-0 
lines. However, we also require that the Q-branch 
lines are not too far below the S-branch lines. In 
this manner, a good approximation for extinction 
is reached in each case. 

4.2. The intrinsic outflow 

The higher K-band extinction in the northern 
lobe by ~ 0.7 mag implies that this lobe is not in- 
trinsically weaker than the southern lobe in molec- 
ular hydrogen emission. The north-south 1-0 S(l) 
ratio was found by Eisloffel et al. (1996) to be 0.54. 
This suggests that the outflow may be much more 
symmetric than appears. Since the northern lobe 



is redshifted, this is also consistent with the lobe 

dynamics within a spherical cloud. 

The mean density can be estimated on assum- 
ing the geometry of the enveloping cloud. Taking a 
uniform spherical cloud of radius 30", correspond- 
ing to 0.11 pc, and an average column to the out- 
flow of ~ 3 X 10^^ cm~^, yields a mean (H nuclei) 
density estimate of 1.3 x 10'''cm^''. 

The extinction also implies that the intrinsic H2 
luminosity from the 1-0 S(l) line alone is 0.33 Lq, 
and the total H2 luminosity will probably lie in the 
range 3-9 Lq, depending on the type and strength 
of shock waves involved (Smith 1995). Notably, 
this total is very close to that estimated for the 
complete far-infrared cooling by Giannini et al. 
(2001) of 2.8 Lq, which includes CO, H2O, O and 
OH emission. 
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5. Ha and CO Modelling 

To simultaneously model the KSPEC, ISO- 
CAM and LWS CO data, we require a calibration 
of the KSPEC data set. We employ the integrated 
1-0 S(l) K-band fluxes of 1.45 x IQ-^^ Wm'^ 
(North) and 2.68 xlO^^'"^ Wm-'^ (South) tabulated 
by Eisloffel et al. (1996) since the associated areas 
correspond quite closely to the effective apertures 
of the ISO observations. 

We have attempted to apply simple planar J- 
shock and C-shock models, without success. We 
have employed an updated shock code described 
in detail by Smith et al. (2003). The code as- 
sumes that a shock is stable and in a steady state, 
the ion number is a conserved quantity and that 
the H2 dissociation rate is given by equilibrium 
conditions. We take an ortho-para ratio for H2 of 
3. Single J-type shock waves predict high excita- 
tion H2 spectra and C-type shocks predict quite 
constant excitation across a wide range of upper 
energy levels. 

Multiple shock waves are required, as demon- 
strated in Fig. 4. Two carefully chosen C-type 
shocks provide a reasonable fit to the full set of 
data. The top panel displays both the ISO H2 
data (stars at low T^ ) and the KSPEC data. The 
lower panel displays the CO rotational fluxes. It 
is not clear, however, how two shocks with almost 
the same parameters can be found in the three 
separate locations since the excitation produced 
in planar shocks is very sensitive to the shock ve- 
locity, field strength and ionization fraction. 

Next, we fit a single paraboloidal C-type bow 
shock. C-type and J-type bow shocks are modeled 
according to a scheme illustrated in Fig. 5. Cylin- 
drical coordinates (z, R, (p) are used, and the mag- 
netic field B, is defined via the density and Alfven 
speed. The spectroscopic results presented here 
arc independent of the direction of the observer, 
a. The molecules axe completely dissociated at the 
bow cap provided the bow is moving faster than 
the appropriate dissociation speed limit. The data 
provides a large number of constraints. For exam- 
ple, the CO fluxes and excitation determine not 
only the CO abundance but also the density. The 
excitation state of the vibrationally excited H2 
also determines the density as well as the shape of 
the bow. The rotationally-excited H2 determines 
the atomic hydrogen fraction. 



The final results are displayed in Fig. 6 for Cep 
E South and Fig. 7 for Cep E North (location Nd 
for the KSPEC data). Note that bow configura- 
tions are strongly supported by two independent 
observations. First, the H2 images display nu- 
merous bow-shaped structures in both lobes (e.g. 
Fig. 1 and Ladd & Hodapp (1997) and, second, 
many locations within these 1-0 S(l) bows pos- 
sess double-peaked line profiles, as predicted by 
bow shock models (Eisloffel 1997) and munerical 
simulations Suttner, Smith, Yorke, & Zinnecker 
(1997). 

The displayed CBOW model fits are further ev- 
idence that these are bow shocks. The ranges of 
critical parameters determined during the mod- 
elling procedure are given in Table 4. Most signif- 
icant are (i) the density, (ii) the molecular fraction, 
(iii) the bow shape and (iv) the CO abundance. 
Less critical parameters to the gas excitation state 
are (i) the bow speed (the excitation is fixed pro- 
vided the bow speed exceeds a minimum value), 
(ii) the Alfven speed (i.e. the magnetic field), (iii) 
the ion fraction (a minimum value is necessary for 
a J-shock to be appropriate) and (iv) the oxygen 
abundance. The latter parameters infiuence the 
location of the molecular emission along the bow 
surface but, provided the location exists and is not 
at the apex, the excitation is not strongly affected. 
Perhaps the most significant result is the atomic 
hydrogen fraction of 0.1 - 0.4. Additionally, the 
derived CO abundances at all three locations lie 
close to 1 X 10^^. 

Finally, we consider the type of J-type bow 
shock that could fit the spectroscopic data. Since 
J-shocks heat the gas impulsively, higher vibra- 
tional levels are relatively well populated. Hence, 
to fit the data, a bow with longer cooled wings 
than a paraboloid is needed to compensate. After 
varying the critical parameters. Fig. 8 shows that 
a bow with shape s = 1.5 where z/L = (l/s)(R/L)* 
fits the data extremely well. 

A major reason that quite blunt shapes are pre- 
dicted here, in comparison to those derived for 
Cepheus A, is that lower densities are involved. A 
lower density is deduced from the relatively low 
CO fluxes in Cepheus E. The lower density then 
necessitates a higher fraction of hydrogen atoms 
to maintain a small divergence from local thermo- 
dynamic equilibrium between the populations of 
the H2 vibrational levels. 
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Table 4: Ranges for the critical parameters which describe the observed shock excitation, as derived by 

comparing models to both the H2 and CO data. Location Nd is taken here. The criterion taken as a good 
fit is a discrepancy of approximately 30% between the majority of data points and the model fluxes. 



Model 


Bow 


Density 


Magnetic 


Alfven 


Ion 


e(C) 


e(0) 


bow 


H2 




speed 


10^ 


field 


speed 


fraction 






shape 


fraction 




km 


cm-^ 


mG 


km 


10-7 


10-4 


10-4 






C-type 


> 38 


0.6-1.2 


0.2-1.3 


1.0-6.0 


4-50 


0.6-1.4 


1.0-5.0 


1.9-2.3 


0.32-0.40 


J-type 


> 17 


2.0-5.0 


0.3-0.6 


1.5-3.7 




0.3-0.6 


1.3-7.0 


1.4-1.6 


0.14-0.26 



6. Further diagnostics 

To distinguish between J and C-type bow mod- 
els, we examine other signatures. Wc present the 
contributions to the cooling, calculated from the 
two best fitting models to Cep E SI, in Table 5. 
The luminosities are calculated by assuming a suf- 
ficient number of bow shocks to explain the H2 
1-0 S(l) luminosity. Indeed, a reasonable number 
of distinct bow shocks are inferred. We find from 
the above model that a single C-type bow of speed 
Vfeow = 120 km s~^ and scale size L = 10^^ cm 
emits 0.016 L© in the H2 1- S(l) line. There- 
fore, 20 such bows in each lobe would suffice to 
provide the observed luminosity. Note that these 
bows would appear larger than L in size, with H2 
being thermally dissociated until the normal com- 
ponent of velocity falls below 45 km s~^. For a 
paraboloid, it is straightforward to show that this 
occurs at R = 2.47 L and z = 3.05 L. 

A single J-type bow with s = 1.5, Vbow = 
60 km s^^ and scale size L = 10^'"' cm emits 
0.02 Lq. Hence, 16 such bows would be needed to 
account for the H2 emission. The apparent bow 
size is, however, larger. With a dissociation speed 
limit of 25 km s~^, we find molecules survive only 
for bow locations at R > 4.76 L and z > 4.36 L. 

The [O I] 63/im line is potentially a good shock 
diagnostic. However, the column density of hy- 
drogen nuclei which provides unit optical depth at 
line centre is just 5 x 10^° cm^^ (HoUenbach & 
McKee 1989). This may explain why the model 
line fiux is overpredicted by factors of a few (see 
Table 5). Given gas columns of 3 x 10^^ cm~^, 
the observed [01] emission may well not have a 
shock origin. Furthermore, we also detected the 
[ClI] (158 iim) line in the LWS spectrum. The mea- 
sured ratio [Ol](63/xm)/[ClI](158yum) is about 1.3 
(1.5) for the north and south lobes, respectively. 



Table 5: Contributions to the intrinsic cooling of 
the flow in Cep E South and the bow shock pre- 
dictions. The total cooling estimates are extracted 
from Giannini et al. (2001). 



Coolant 


Observed 


C-Bow 


J-Bow 




Lq 


Lq 


Lq 


H2 1-0 S(l) 


0.33 


0.33 


0.33 


H2 dissociative 




0.11 


0.92 


H2 rot+vibrat. 




6.77 


8.32 


CO rote 


>0.51 


4.76 


56.56 


CO rote (J > 10) 


0.51 


1.26 


1.38 


CO vibe 




0.01 


0.03 


H2O 


>0.67 


1.53 


3.71 


H2O (FIR, T>800K) 


0.67 


0.54 


1.39 


OH 


0.17 


0.01 


0.01 


Ol 63/im 


>0.21 


1.72 


21.36 


CI 370/im 







0.17 



In shocks, the [Cll](158/um) line is usually several 
orders of magnitute fainter than the [Ol] (63 /im) 
HoUenbach & McKee (1989). Hence these lines 
may have a PDR origin, as already discussed by 
Moro-Martfn et al. (2001). The observed OH lu- 
minosity is, according to Giannini et al. (2001), to 
be treated with caution since it may well be due 
to continuum pumping. 

We conclude from Table 5 that the long wings 
of the J-type bow produce extremely high fluxes 
in the lines from the cooler molecular gas. This 
is inconsistent with all expectations since the to- 
tal cooling then exceeding the bolometric lumi- 
nosity of the protostar. Furthermore, in a shock- 
dissipative momentum-driven outflow, we expect 
the mechanical power of the flow to be very close 
to the total cooling. We estimate here the mechan- 
ical power through CO J=2-l observations by first 
noting that most of the emission is from low-speed 
CO but most of the kinetic energy lies in the CO 
moving with radial speeds of 15-25 km (Smith, 
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Suttncr, & Yorkc 1997) (a result of the very flat 
line profile found for Cepheus E). Hence, given a 
blue lobe mass of 0.16 M© Ladd & Hodapp (1997), 
we derive a mechanical energy of ^ 5 x 10*''' erg 
and a mechanical power ~ 20 Lq. Although this 
estimate is subject to considerable uncertainty, it 
is consistent with the C-type bow model. 

A further diagnostic is provided by the inte- 
grated intensity of the CO (3 2) emission shown 
in Fig. 2. The C-bow model predicts a C0(3 
2)/CO(18 - 17) ratio of 6.0 in contrast to the J- 
bow model (276) for the southern outflow lobe. 
The measured ratio is 4 (north) and 2 (south), 
much nearer to the predictions of the C-type bow 
model. 

The columns in the 0-0 S(3) and 0-0 S(5) ap- 
pear lower than predicted in the southern lobe. 
This suggest that the ortho-to-para ratio of H2 
may be under three, the LTE value for high tem- 
peratures. This can occur in C-shocks since the 
gas is gradually heated from a low pre-shock tem- 
perature. Therefore, the gas maintains the pre- 
shock ortho-para ratio. After some heating, how- 
ever, atomic hydrogen can induce ortho-para con- 
version (Smith, Davis, & Lioure 1997). 

Note that the distribution of shocks generated 
by a supersonic turbulent velocity field can pro- 
duce an excitation that mimics that from a bow 
shock surface (Eisloffel, Smith, & Davis 2000). 
It can be shown that a power law number dis- 
tribution of shock speeds (with power law index 
a = — 2s/(s — 1) where N oc f~") is required. This 
also assumes that there is no systematic magnetic 
field effects although, for bow shocks with moder- 
ate Alfvcn speeds, we find that the magnetic field 
direction has an insignificant influence on the ex- 
citation. 

7. Conclusions 

We have analysed the molecular outflow from 
the Cepheus E-MM source over a broad wave- 
length range, from the near infrared into the 
sub-millimeter regime. We combined KSPEC 
NIR data, ISO mid- and far-infrared spectra and 
JCMT sub-mm observations. We have concen- 
trated on the properties of the vibrationally and 
rotationally excited H2 and CO since sufficient 
data points exist to permit a simultaneously in- 
terpretation in terms of shock models using J- 



and C-type physics. We investigated planar and 
curved shocks. 

Cepheus E is quite exceptional in its strong ra- 
diative power (above 10 L©) and short dynamical 
age (under 700 yr given an advance speed of over 
100 km s~^). Our main results are as follows: 

1. Extinction is most accurately derived by sat- 
isfying multiple constraints on the complete 
NIR data set rather than from specific line 
ratios. The southern blueshifted lobe has 
a mean K-band extinction of 1.4 mag and 
the northern redshifted lobe has extinction 
in the range 2.1 -2.4 mag. This is consis- 
tent with the extinction being internal to a 
spherical cloud, with a density of 10^ cm~^ 
and radius 3 x 10^^ cm~^. Remarkably, the 
extinction-corrected H2 columns, including 
the ISO data, demonstrate a state close to 
local thermodynamic equilibrium. 

2. We interpret all the measured H2 and CO 
line fluxes simultaneously in terms of shock 
models. The best fitting models are obtained 
using shock distributions in the form of bow 
shocks. 

3. C-type physics is strongly favoured mainly 
because J-type physics predicts extremely 
strong emission from the low-excitation 
flanks of a bow, which is not observed. 
Moreover, for the J-shock model, the outflow 
power and momentum become implausibly 
large in comparison to the values found for 
the protostar and the bipolar outflow, re- 
spectively. 

4. High resolution H2 1 S(l) imaging has 
shown that the lobes consist of numerous 

emission knots of size 1 3 x 10^^ cm (Ladd 

& Hodapp 1997). We find that of order 20 
such bow shocks, close to paraboloidal in 
shape, are required by the model to explain 
the intergrated emission for each lobe. 

5. No significant differences in the density or 
abundances are found between positions. 
We also find no significant differences in the 
H2 excitation when analysing the integrated 
spectra of a whole outflow lobe or a 6" x 6" 
pixel-sized part. This is in agreement with 
Eisloffel et al. (1996) and suggests that the 
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large bows arc built up of smaller unresolved 
bow shocks, generated by flow instabilities. 

6. The pre-shock medium is not fully molecu- 
lar in these models. We have found a mean 
atomic fraction of n{H) / {n{H) + 2n(i?2)) = 
0.4. With bow speeds of 120 km s~^, exten- 
sive bow apices are predicted within which 
molecules are completely destroyed. The 
reformation time is of order 10^^/nc(H), 
where nc(H) is the atomic density in the 
compressed layers. Hence, a reformation 
time of ~ 3000 yr may be achieved for 
nc(H) = 10^ cm-3. 

7. The total mass set into motion by the out- 
flow is low, 0.25 M0 (Ladd & Hodapp 1997). 
We estimate a total cloud mass of 5Mq 
given a density of 10^ cm~^ and radius 
2.2 X 10^^ cm. With an outflow half-opening 
angle of /?, a fraction fc = 1-cos/? of the 
cloud would be disturbed by two lobes. Tak- 
ing f3 = 20°, yields fc = 0.06, and thus pro- 
vides a consistent basic model. 

High outflow powers have been uncovered from 
many other Class protostars. Mechanical lumi- 
nosity exceeds 50% of the bolometric luminosity 
for all of the Perseus Class sources (Barsony ct 
al 1998). On the other hand, a few Class pro- 
tostars such as L1527 and B335 possess relatively 
weak outflows as niciasured by total far-infrared 
luminosities (Giannini et al. 2001) and CO mo- 
mentum flow rates (Bontemps et al 1996). In this 
respect, the Cephcus E outflow could represent a 
powerful but abrupt evolutionary phase, about to 
be brought to a halt as underlying jets exit a com- 
pact 0.1 pc cloud. 

Jochen Eisloffel and Dirk Froebrich received fl- 
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forschung grant 50 OR 9904 9. The ISO Spectral 
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Fig. 4. — Planar shock models for CepE South. 
In addition to the symbols described in Fig 3, 
the ISOCAM data are represented by stars, (a) 
A hot C-shock component with 32kms~^ pro- 
vides a fit to the vibrationally-excitcd columns 
but does not contribute to the CO rotational 
lines (other critical parameters are a density of 
10^ cm~^, an Alfvcn speed 1.5 km s^^, ion frac- 
tion X 10^^, a transverse magnetic field and 
n(H2)/(n(H)-h2n(H2)) = 0.3). The maximum 
temperature is 3160 K. (b) A cool C-shock with 
speed llkms"^ provides a fit to the rotationally- 
excited H2 and CO but does not contribute to the 
vibrational lines (other critical parameters are as 
above except n(H2)/(n(H)-h2n(H2)) = 0.49998). 
The maximum temperature in the shock is 950 K. 
Also modelled are abundances of e(0) = 4 10~^ 
and e(C) = 1.5 10"^. 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.0. 
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Fig. 5. — The geometrical parameters associated 
with a fast-moving bow shock, as originally em- 
ployed by Smith and Brand (1990). 



Fig. 6. A C-typc bow shock model for 
CepE South 1. A C-type bow model with 
speed 120kms~^, a density of 10^ cm~^, an 
Alfven speed 1.5 km s~^, ion fraction x 10~^ 
and n(H2)/(n(H)+2n(H2)) = 0.35, abundances of 
e(0) = 3 10-^ and e(C) = 1.4 10"^ (initially 
in CO) and a magnetic field aligned with the bow 
axis. The H2 model values are displayed for the 
ground (solid line ) , first (dashed) and second (dot- 
dash) vibrational levels. 
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Fig. 7. — A C-type bow shock model for CepE 

North (location Nd). The model bow speed is 
again 120kms~^ and critical parameters are as in 
Fig. 6, except for the initial abundance e(CO) = 
1.2 X 10^^. Displayed lines on the CDR diagram 
represent individual vibrational levels (see Fig. 6). 
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Fig. 8. A J-type bow shock model for 
CepE North, position Nd, with shape parame- 
ter s = 1.5. A speed of 60kms~^, a density 
of 3 10'* cm~^, an Alfvcn speed 2.5 km s"-'^, 
n(H2)/(n(H)+2n(H2)) = 0.2, and abundances of 
e(0) = 4 10-^ and e(C) = 5 10"^ (initially 
in CO) and a magnetic field aligned with the bow 
axis were taken. The H2 model values are dis- 
played for the ground (solid line ), first (dashed) 
and second (dot-dash) vibrational levels. 
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